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Abstract 



O . 

O" We prove a removable singularities theorem for stationary Einstein 

£h . equations, with useful implications for constructions of stationary so- 

lutions using soliton methods. 



(N 1 Introduction 



The soliton technique has proved to be a very effective tool for construct- 
ing stationary black holes in five dimensions, see e.g. [Il[6j[7J[T6j[I7]. The 



method is used to construct singular solutions of harmonic-map type equa- 
tions. One then needs to make sure that the singularity structure of the 
resulting harmonic map is compatible with a smooth geometry of the asso- 
ciated space-time. 

rN | Proceeding in this way, in their ingenious construction of Black Sat- 

ed ' urns [3], Elvang and Figueras introduce a singular point 

«i := (p = 0, z = Oi) 

on the boundary {p = 0} of the Weyl coordinates domain {p > 0, z G R}, 
and fine-tune certain constants to ensure that the metric functions remain 
uniformly bounded near a±. Now, the resulting metric ends up being a 
rational function of 



R\ ■= VP 2 + (z-ai 
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This leads to a potential problem because Ri is not differentiable at a\, and 
therefore the differentiability of the metric functions at «i is not apparent. 

We will refer to such points as ghost points, as their occurrence does not 
seem to be related in any obvious way to desirable geometric properties of 
the resulting space-time such as end points of horizons or fixed points of the 
isometry action, compare [TUHT2"] . 

Closer inspection [2] of the Black Saturn metric near a\ shows that, with 
the choice of free constants that makes the metric functions bounded, all 
the metric functions can be rewritten as rational functions of R\; since this 
last function is smooth, smoothness of the metric near ot\ becomes obvious. 
The calculations required to establish this fact turn out to be rather heavy, 
requiring quite a bit of effort to coerce Mathematica to produce the re- 
sult. We emphasize that the result is non-trivial and requires non-obvious 
factorisations and cancellations of odd-order polynomials in R\. 

A similar trick of introducing ghost points has been used in other related 
constructions [5|l8llT3 |[T8llT9] . The question then arises, whether one needs to 
redo the calculations of [2] case by case, or there is a general mechanism which 
guarantees that ghost points are smooth points for the resulting metric. 

The object of this note is to show that smoothness of the metric at such 
points is a consequence of the stationary Einstein equations with matter 
fields, without the need to assume more Killing vectors. This can be roughly 
stated as follows, see Theorem 12.11 for a precise version: 

Theorem 1.1 Singularities of Lipschitz continuous stationary Einstein met- 
rics located on timelike submanifolds of codimension m > 2 are removable. 

We discuss in somewhat more detail in Section [3] how this theorem takes 
care of the ghost point problem. 

2 Stationary Einstein equations 

We consider a time- independent metric in a space-time of dimension n + 1. 
Since the problem we address is purely local, we assume that the space-time 
metric functions q^ are given in local spatial coordinates ranging over a ball 
B{R) C M. n , n > 2, of radius R centred at the origin, and the prospective 
singularities lie along a smooth submanifold 



of B(R) of codimension 2 < m < n, with either <9£ = or dT, C dB(R). We 
set 

B*(R) :=B{R)\H . 

In adapted coordinates the metric can be written as 

= -v 2 (dt + W) 2 + gadyW , (2-1) 

=8 ~xr 

where d t is (stationary) Killing, i.e. 

dtV = d t 6 = d t g = 0. (2.2) 

We allow matter fields ip = (<p A ) with energy-momentum tensor that de- 
pends upon g, dg, <p and dip. For simplicity we assume that the <p A, s trans- 
form as scalars or tensors under coordinate changes, and that the stationary 
matter field equations constitute a tensorial system of the form 

A g cp = F(g, dg, cp, dip) in B*(R) , (2.3) 

though a wider class of more general elliptic systems can be easily incor- 
porated in our analysis. We note that (linear) electromagnetic fields, for 
example, satisfy this assumption in Lorenz gauge. 

The Einstein equations with (possibly zero) cosmological constant A for 
a metric satisfying fl2.ip - fl2.2p read (see, e.g., [3] or pQ) 

rl 



div 9 (VA) = 2V T - ( ^±±A - %gl ) V 2 9 ) (2.4) 



' VA q V = -l\\\ 2 + Too-(^A-^^)v 

y 41 Iff uu l n-l re— 1 I 

div 9 (V/A) = 2v[t - (giA - t ±g 1 )v 2 6 
^ Ric(<?) - ^Hess^V = ^A o A + T, + (Sl A " % i T 1 )^ > 
where T g := T i jdx 1 dx :) , T = T 0i cfe J , and 

Ay = -y 2 ^^ - a,^) , (a o \) tJ = \ k \ k] . 

Altogether the Einstein-matter field equations, which are supposed to hold 
in B*{R), can therefore be written as 

A g ip = F(g,dg,<p,d<p) , 

A g V = F 1 (g,dg)+f v , 

dw 9 (d6) - d(dw g (6)) = F 2 (g, dg) +T e , [ ^ 



Ric(g) - V-tResSgV = F 3 (g, dg) + T t 
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for some (explicitly computable) Fi, F 2 and F 3 which are polynomials in g, 
_1 , dg and quadratic in 8q and some TV, Tg and T g which arise from T^. 

We have: 

Theorem 2.1 Under the conditions above, suppose that (g, ip) G C°' a (.B(.R))n 
C 2 {B*(R)) andV > ^n £(i?). Assume further that 

1. either a = 1, 

2. or "^—^ < a < 1 and t/iere exists a constant C, possibly depending upon 
(g, p), such that 

Idsl + I^I^CdistR^E)"- 1 , (2.6) 

and 

\T^\ + \F\<C(l + \dg\ 2 + \d<p\ 2 ). 

Then 

B ,<peC u (B(R/2)). 

The proof will use the following simple lemma, whose proof is deferred 
until after the proof of Theorem 12.11 

Lemma 2.2 Let Q be an open subset ofM. n and £ C fi be a smooth subman- 
ifold of codimension m > 2 which either has no boundary or has boundary 
contained in dQ. Assume that u G W l,2 (Q) satisfies 

d i (a ij d j u) = d i g i + f inQ\E 

in the sense of distributions for some a v -> G L°°(Q), f G L 1 (fi) and g l G 
L 2 (Q). Then u satisfies the above equation in Q in the sense of distributions. 

Proof of Theorem 12. 11 By hypothesis there exist coordinates y % in which 
the metric coefficients and the fields satisfy 

V,9 ugij ,p e C°< a (B(R)) nC 2 (B*(R)) . (2.7) 

Standard arguments (compare [15]) show that the metric is in fact smooth 
away from E: 

V,e t ,g l3 ,p eC^ a {B{R))r\C°°{B\R)) . (2.8) 



Since the problem is local, it suffices to establish the desired regularity 
in a small ball B(e) centered at a point on E n B(R/2), which is assumed 
to be the origin. By a linear change of coordinates we can without loss of 
generality assume that gij(0) = 5?. Furthermore, we can also assume that 

E n B(10e) = {y = (y\ 0) : y' e tt(E)} (2.9) 

for some open set tt(E) C S(10e) n M n_m 

Most of our arguments rely on elliptic estimates. Of the four equations 
in (12. 5p . the last two are not manifestly elliptic. As is well known, this 
issue can be cured by passing to harmonic space-coordinates and using an 
appropriately chosen time function. 

Suppose, first, that ^^ < a < 1. 

For e > let x l be solutions of the problem 

A g x i = 0, x i \ S (e) = y i . 

where S(e) := clB(e) is a ^-coordinate sphere of radius e. 

By [91 Theorem 8.34], and elliptic regularity away from the origin, we 
have x 1 G C l ' a (B(e)) n C°°(B(e) \ E). If we write 

x l = y l + r , 

then f l solves the divergence-type equation 

dii^/d^g^djf) = -d^y/MTgg 11 ) = -d^^d^g^ - 8f) . 
Let 



f!(x) := f\ex) , tf(x) := y/tetg{cx)g u {ex) , $?{x) := tf(x) - tf (0) . 

Then 

d i (^d j fi) = -ed i ^ in 5(e). 

Applying P, Theorem 8.33] to the equation satisfied by f[ on a ball of radius 
two we obtain (note that the first term |w|o there can be discarded by the 
usual argument that exploits injectivity of the Laplace equation) 

\\f e \\c^(B(i)) < Ce ||ffico.-(B(i)) < Ce 1+a . (2.10) 

It thus follows that 



dx 
~dy~i 



— = S^ + 0{e a )mB{e). (2.11) 



The implicit function theorem shows the x*'s can be used as a coordinate 
system near y % = for all e small enough. In what follows we choose some 
such value of e. 

As for the choice of a spacelike slice, we will now show that we can assume 
without loss of generality that 

dw g (9) = in B(e) for sufficiently small e . (2.12) 

To use the freedom of defining t (equivalently, to fix the gauge freedom of 9), 
we make a coordinate change of the form t = t + h(y). The metric q then 
takes the form 

= -V(di+ (9i - hi) dy 1 ) 2 + 9ij dy l dy j . 

To obtain fl2~T2l) . we pick h e W 1 ' 2 (B{e)) to be a solution to 

A g h = div g (9) in B(e) , and h\s( e ) = . 
By [HJ Theorems 8.12 and 8.34], we have 

h e C l ' a {B{e)) n W 2 ' 2 (B(e)) n C°°(B(e) \ E) . 
Recalling (12.61) and rewriting the above equation for h as 

g ke d k d e h = -d t hd k (^g~g M ) + -L=d k {y/dtid Q k %) , 

\/detg 

we can apply (HI Theorem 9.11 and Lemma 9.16] to get, 



d 2 h e L q (B(e)) for any 1 < q < 



in 



1 — a 



here and in what follows the norm might depend upon e, but this is irrelevant 
since a small epsilon has been now fixed. We have thus achieved (I2.12p with 
a penalty that the derivatives 89 no longer satisfy a pointwise estimate given 
by ( 12. 61) but the weaker estimate 

IT) 

89 e L q (B(e)) for any 1< q < . (2.13) 

1 — a 

The above suffices for our purposes. We emphasize that V and g remain 
unchanged under the redefinition of 9, equivalently of time, as above. 
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To prepare for our passing to the coordinates x l , we need some bound for 
the Hessian of /\ For y e W 1 , we will write y = (y',y") where y' G M n_rn 
and y" e W 71 . In view of (J2.9P we have 



d(y) := dist R n( 2/ , E) = |j/"| for 3/ e £(5e) \ £ . (2.14) 



Define 



Q' s ,s = {y ■ s/A < \y"\ < 5s/4, \y\ < S} , 
Q s ,s = {y : s/2 < \y"\ < s, \y\ < S} . 

Then, by ( 12. 6p . /f satisfies 

\^ s j d.dj f s \ < Cs a+1 in Q[^ le . 

Thus, by Theorem 9.11] and ( 12~T0|) . 

||5 2 /ilU,(Q 1>s _ l£/2 ) < Cs a+1 [H"" m (S s )] V,? for any 1< q < oo , 

where '}i n ~ rn denotes the (n — m)-dimensional Hausdorff measure: More pre- 
cisely, we first apply [91 Theorem 9.11] to cubes of unit size and fl — L(f^) 
with L(fg) being the linearization of /* at the center of those cubes, and 
then sum the acquired estimates over a collection of non-overlapping cubes 
covering the desired region. Because of the simple geometry of S s (compare 
( 12. 9p ). the number of cubes in each such collection is proportional to s m ~ n , 
which is itself proportional to the Hausdorff dimension above. Scaling back, 
it follows that 

\\d 2 f\W{Q s , e/2 ) <Cs a - 1+m /« for any 1< q < oo . 

Now if we pick q such that a — 1 + m/q > 0, we can sum the above over 
dyadic rings in the tranverse direction to get 

771 

\\d 2 f\\ L HB(e/2)) < C for any 1< q < . (2.15) 

1 — a 

We pass now to the coordinates x % = y % + /*, and still use the symbol g, 
g, 9 and <p for the space-time metric, the spatial metric, the shift one-form 
and the matter fields in the new coordinates. Shifting the x J 's by a constant 



vector if necessary, we can assume that x l (0) = 0. Furthermore, one has in 
the new coordinates 

F,e i)ft ^eC°' a (B(6))nr(B*(e)). 

Estimate (123511 shows that x i G W 2 ' q (B(e/2)), and by (J2SJ, fl233|) . the 
chain rule and the transformation law for tensors one deduces that 

Tfl 

\d x g\ + \d x <p\ G L q (B(e/2)) for any 1< q < - . (2.16) 

1 — a 

In the coordinates x % the Einstein-field equations ( 12. 5 j) can be rewritten 
in the following form 



A fl y = F 1 (fl,9g)+2V, 

A/ i = F (i) (g,9 S )-(f 9 ) i , 

Ap^j - 9j(9j log y) = F(i)(j)(0, dg) + (f g )ij , 



(2.17) 



where we have used ( 12.121) . Using ( 12.161) together with the given growth rate 
of F and T, one sees that (V,9,g,<p) G W l,2 (B(e/2)) while the right side of 
(12371) belongs to L p (B(e/2)) for any p < ^fz^y- Also, by Lemma E2J ( |237j) 
is satisfied across S in the sense of distribution. It is useful to write the last 
equation in (12.1 T[) as 

&g9ii = d i( d J lo S V) + F 3 ( S , d B ) + (f g )ij . (2.18) 

To proceed, we distinguish two cases according to whether a > 1 — 7p- 
or a < 1 — y~. In the former case, we apply [HJ Theorem 5.5.3(b)] to 
the first three equations of (1237J) to assert that (<p,V,6) G C 1 ' a (B(e/3)) for 
some a > 0. In particular, <91ogU G C ,0,cr (i?(e/3)). Applying p3| Theorem 
5.5.3(b)] again to (123HJ) . we get G C7 1,<T (5(e/4)). 

In the latter case, we use [HI Theorem 5.5.3(a)]. Applying this result 
to the first three equations in ( I2.17P and then to (I2.18P as in the previous 
paragraph we get (</?, V, 9,g) G W 1,9 (S(e/4)) for any 1 < q < 2(1 _^_m . In 

other words, in B{e/A), (12.161) is improved with a replaced by a + (a — ! ^ iri ) • 
Repeating this process for a finite number of time, we arrive at a situation 
when the argument in the previous paragraph applies. 

In any event, one obtains (g,<p) G C 1,o "(.B(e/100)). A standard bootstrap 
argument based on Schauder estimates proves smoothness; analyticity readily 
follows. 

8 



When a = 1 we replace a by any number in (1 — ^,1) and arrive to (12.161) 
as before. Since T^ is a tensor, it gives a bounded contribution to (12.1 7ft . 
leading to a metric with improved regularity as before. Similarly F is a tensor 
giving a bounded contribution to (12. 3p . and we obtain (g, <p) e C 1,a (B(e/2)) 
by the same method as above. The result follows. □ 

To finish this section, we provide the 

Proof of Lemma 12.2b Let £ e C™(Q), we need to show that 

f a ij d i ud£dx = I ' (g* d& + f £) dx . (2.19) 

Ju Ju 

Let rj be a smooth cut-off function on R such that r/(t) =0 for t < 1 and 
r](t) = 1 for t > 2. For 5 sufficiently small, define 

77 ( ? ) for ?ti > 3 , 



aWH ^4s^)fo"« = 2; 



log(-log(5) 

By hypothesis we have 



/ a lJ diU dj(£ Q 6 ) dx = (g l dt(£ Qs) + f £ Qs) dx . 

Jn Jn 

(I2.19P can be reached by passing 5 — )■ using Lebesgue's dominated con- 
vergence theorem, Cauchy-Schwarz's inequality and the explicit form of g$. 
Note that when m = 2, we need to use 

dt < 00 . 



tilogtf 
We omit the details. □ 



3 Ghost points 

We show how Theorem 12.11 applies to stationary solutions obtained by intro- 
ducing ghost points in the solitonic solution-generating technique. Here the 
space-time metric that one wishes to construct is invariant under an abelian 
isometry group M x T n-2 , where the R factor represents ^-translations. The 

9 



metric depends only upon two coordinates (p, z), which can be thought of 
as cylindrical coordinates on M. 3 . In this construction, if the ghost point is 
placed at (p — 0,z — 0), one obtains a solution of the vacuum Einstein 
equations g^ defined on 

{t6l, (x, y, z) e B{8) C M 3 } x T n " 3 , 

where 

x = p cos ip , y = p sin ip . 

Note that one S 1 factor from T™ -2 has been interpreted as a rotation around 
the z-axis of 1R 3 . Furthermore, there exists a neighborhood of the origin in 
which the metric functions are analytic functions of (p 2 , z, d), where d is the 
Euclidean distance to the origin in IR 3 : 



d := V p 2 + z 2 . 

So the singular set 

Y 1 = { x = y = z = 0}x T n ~ 3 

has dimension n — 3 within each slice t = const, hence codimension three. 
To apply Theorem 12. II we need to verify that, reducing S if necessary, 

there exists e > such that detjf- 7 > e and q u < —£• (3.1) 

Note that the function d is Lipschitz-continuous, but not differentiable. This 
implies that the metric functions are in C°°(B*(R)) fl C Q,1 (B(R)). Theo- 
rem 12.11 with a = 1 shows then that the metric functions are real-analytic in 
a whole neighborhood of the origin of M. n , as desired. 

In the case of the Black Saturn metric we have n — 4, and the space- 
dimension of the singular set is one. To verify that this metric is analytic 
near its ghost point a.\ = (p = 0,z = ax), one needs to verify (13. ip . A 
direct verification in the coordinate system used in [1] fails, because at this 
point Q u becomes null for all Black Saturn metrics (indeed, a\ always lies on 
the ergosurface for those metrics). This can be bypassed by checking that 
the limit of the metric at a.\ is Lorentzian, and that the determinant of the 
matrix of scalar products of all Killing vectors there has a strictly negative 
value. This guarantees that some linear combination of Killing vectors is 
timelike at a±, and our theorem applies. 

As another application, our analysis reduces the question of regularity 
of the metrics of [51IS1H31HH1HH] near their ghost points to showing that the 

10 



components of the metric tensor in coordinates (x, y, z) have a finite limit 
at the ghost points, and verifying (13.1 j) there (after perhaps replacing d t by 
a different, timelike Killing vector if necessary). We note that this might 
require tedious symbolic algebra calculations, and our experience with the 
Black Saturn metrics suggests that the checking of the timelike character of 
the orbit of the isometry group through the ghost point might be non-trivial. 
In any case we have not attempted to carry this out. 
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